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A* Input high arxJ low energy- 
Imases. Shrink them to a 
smaller matrix size. 
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B. Compute gradient mask using 



.204 
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LUT 

Containing 
mapping of 
w c 's and w- 
ranges for 
high and 
low energy 
kVp pairs 



C. Localize mosl 
eliminate ribs usin 
constraint on the 1 
as intensity thresh 


erf me bones and 
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— H. Refine We** =Wb 
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Figure 2 
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D. At each point on me 
gradient mask, vary w over the 
range and obtain w corresponding 
to the minimum gradient 
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£. Construct the w-histogram 
of aK the points in the mask 
and find the peak to obtain the 
maximum Ifcellhcod estimate Wml 
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F. If peak is at the extreme of 
range, 

Else w»=Wml 
wife s a(wO + b 
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I. Use w» and v*, in Equation 1 on 
the full resolution tow and high 
energy images fo obtain soft tissue 
and bone Images 
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Compute decomposed image I(x,y) utilizing the largest parameter in 
the cancellation range w-range (e.g. w M emory), the high energy level 
image 116 and the low energy level image 1 18 in Equation 1 
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Compute gradients (e.g. gradientSMEMORY) utilizing decomposed 
image I(x,y) 



, — 302 



Save gradientSMEMORY and w M emory in memory 
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Compute the next cancellation range w-range by varying cancellation . 
range w-range by .01 (e.g. w S econd = wmemory - .01) 



NO 
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Is wsecond < w-range? 





Stop 
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Compute decomposed image I(x,y) utilizing w S econd, the high energy 
level image 116 and the low energy level image 1 18 in Equation 1 



, 312 



Compute gradients (e.g. gradientSsECONo) utilizing decomposed 
image I(x,y) 



Identify the first pixel (x,y) location to be compared, 
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Is gradientSMEMORY < gradientssEcoND 
at pixel (x,y)? 



YES 



NO 



Save gradientssEcoND anci Wsecond in 
memory. Reassign parameter values 
of gradientssEcoND to gradients M EMORY 
and wsecond to Wmemory 



Is pixel (x,y) the last pixel location? 



NO 



/ — 



324 YES 



Continue to next pixel (x,y) location. 
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Compute the next cancellation range w-range by decreasing wsecond 
by .01 



.326 



Figure 3 
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Figure .7 
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